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Abstract. We present the photometric sample of a faint 
galaxy survey carried out in the southern hemisphere, 
using CCDs on the 3.60m and NTT-3.5m telescopes at 
La Silla (ESO). The survey area is a continuous strip of 
0.2° X 1.53° located at high galactic latitude ~ —83°) 
in the Sculptor constellation. The photometric survey pro¬ 
vides total magnitudes in the bands B, V (Johnson) and R 
(Cousins) to limiting magnitudes of 24.5, 24.0, 23.5 respec¬ 
tively. To these limits, the catalog contains about 9500, 
12150, 13000 galaxies in B, V, R bands respectively and 
is the first large digital multi-colour photometric catalog 
at this depth. This photometric survey also provides the 
entry catalog for a fully-sampled redshift survey of ^ 700 
galaxies with R < 20.5 (Bellanger et al. 1995a). 

In this paper, we describe the photometric observations 
and the steps used in the data reduction. The analysis of 
objects and the star-galaxy separation with a neural net¬ 
work are performed using SExtractor, a new photometric 
software developed by E. Bertin (1996). By application of 
SExtractor to simulated frames and comparison of multi¬ 
ple measurements, we estimate that the photometric ac¬ 
curacy of our catalog is ^ 0.05"* * for R < 22. Then, we use 
a method to obtain a homogeneous photometric scale over 
the whole survey using the overlapping regions of neigh¬ 
bouring CCDs. The differential galaxy number counts in 
B, V, R are in good agreement with previously published 
CCD studies and confirm the evidence for significant evo¬ 
lution at faint magnitudes as compared to a standard non 
evolving model (by factors 3.6, 2.6, 2.1). The galaxy colour 
distributions B—R, B—V of our sample show a blueing 
trend of ~ 0.5”* between 21 < R < 23.5 in contrast to the 
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V—R colour distribution where no significant evolution is 
observed. 


1. INTRODUCTION 

During the past decade, the photometric and spectro¬ 
scopic surveys have allowed to improve our knowledge 
of galaxy formation and evolution. Since the 1980’s, the 
new technology based on CCD detectors has improved 
the photometric efficiency (sensitivity, linearity, high dy¬ 
namic range,...) as compared to photographic plates. Al¬ 
though the photographic plates are un-rivaled in their 
ability to cover large areas of the sky, the large improve¬ 
ments brought by digital surveys, both in the photometric 
accuracy and in the faint limiting magnitude which can 
be reached, have allowed important new insight into the 
properties and the evolution of the different galaxy pop¬ 
ulations. Several deep digital surveys (Tyson 1988, Lilly 
et al. 1991, Metcalfe 1991, 1995, Driver et al. 1994, Small 
et al. 1995) were performed in different regions of the sky 
(with a typical field size < 0.012 deg^) at faint (B< 25) 
or very faint magnitudes (B< 27.5). First, these surveys 
show that the number-counts of galaxies at B ~ 22 are 
in excess with respect to a non-evolving model and the 
disagreement increases with the apparent magnitude. Sec¬ 
ond, the galaxy colour distributions become significantly 
bluer at fainter magnitudes. 

Several models have been elaborated based on the cos¬ 
mological parameters (Hg, qo) and the luminosity func¬ 
tion parameters (L*,<I>*,q;) allowing a good fit to the 
photometric data (see Koo Kron 1992 for a review). 
Models with pure luminosity evolution where L* evolves 
with look-back time (Bruzual 1983, Guiderdoni & Rocca- 
Volmerange 1990, Yoshii & Takahara 1988) predict that 
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the tail of the redshift distribution of very faint galaxies 
should be extended towards high redshifts. This is not ob¬ 
served in the recent redshift surveys to B < 24 (Colless et 
al. 1990,1993, Cowie et al. 1991, Lilly et al. 1991, Tresse et 
al. 1994 and Glazebrook et al. 1995a), which are in good 
agreement with the redshift distribution expected for a 
non-evolving model. Models with number-density evolu¬ 
tion (Rocca-Volmerange & Guiderdoni 1990, Broadhurst 
et al. 1992) are based on a population of dwarf galaxies 
at z~ 0.3 which would have merged into brighter galaxies 
by z~0. However, these models are difficult to reconcile 
with both the recent observations of weak clustering in 
the correlation function of faint galaxies (Efstathiou et al. 
1991 and Roche et al. 1993) and with the physical mech¬ 
anisms for merging (Ostriker 1990, Dalcanton 1993). The 
most recent models use a new estimation of the slope a 
of the local luminosity function (Koo et al. 1993, Driver 
et al. 1994) which is assumed to increase from a ~ —1.1 
to a ~ —1.8 by the presence of a large number of dwarf 
galaxies (dI,dE). This model is supported by the recent 
observations of the Medium Deep Survey with the HST 
to 1=22 (Glazebrook et al. 1994b), where the counts of 
morphologically normal galaxies are well htted by a non¬ 
evolving model and where a large excess of Irregular and 
Peculiar galaxies is detected which could contribute to the 
excess of blue galaxies. 

Eaint photometric and spectroscopic surveys also pro¬ 
vide maps of the distribution of galaxies in three dimen¬ 
sions. The nearby surveys show that galaxies are dis¬ 
tributed within sharp walls delineating voids with diame¬ 
ters between 20 and 50/i“^ Mpc (de Lapparent et al. 1986, 
Geller & Huchra 1989, da Costa et al. 1994) (where Hq 
= 100 h km.s~^.Mpc~^). Very deep pencil beam surveys 
were obtained in particular directions of the sky (Broad¬ 
hurst et al. 1988, Colless et al. 1990, 1993, Cowie et al. 
1991, Lilly et al. 1991, Tresse et al. 1994 and Glazebrook 
et al. 1995a), and in some of these an apparent periodic¬ 
ity on scales of ~ 128h~^ Mpc has been detected (Broad¬ 
hurst et al. 1990). At these depths (z ~ 0.3—0.5), the large 
amount of time required to obtain the redshift distribu¬ 
tion for a complete magnitude-limited sample constrains 
observations to narrow solid angles. Although the existing 
deep pencil-beam probes are adequate for establishing the 
evolutionary history of galaxies, biases caused by sparse 
sampling may affect the data when used to study the large 
scale structures (de Lapparent et al. 1991, Ramella et al. 
1992). 

In this context, a deep redshift survey near the south¬ 
ern galactic pole was started with the main goal to charac¬ 
terize the large-scale structure at large distances (de Lap- 
parent et al. 1993). The spectroscopic survey covers a con¬ 
tinuous solid angle of 0.3 deg^ and contains ~ 700 galaxies 
with R < 20.5 (i.e. B < 22). The median redshift is at z 
~ 0.3. The entry photometric catalogue for the redshift 
survey was obtained by observing in the B, V, R photo¬ 
metric bands up to 24.5, 24.0, 23.5 respectively in a longer 


strip of 0.4 deg^ (this area was not fully covered by spec¬ 
troscopic observations). The photometric data provide an 
adequate sample for measuring with high conhdence level 
the galaxy number counts and the distribution of galaxy 
colours. The description of the spectroscopic sample of the 
ESO-Sculptor survey is given in Bellanger et al. (1995a) 
and the Hrst results about the properties of the large-scale 
structure are shown in Bellanger & de Lapparent (1995b). 

Here we describe in detail the procedures used in the 
reduction and analysis of the photometric sample. The 
paper is organized as follows. In section 2, we describe 
the photometric observations. Sections 3 and 4 outline 
the data reduction procedures and analyses. In sections 
5, 6 and 7 we discuss the transformation into astronomi¬ 
cal coordinates, the photometric calibration of fields and 
the magnitude transformations into the Johnson-Cousins 
standard system. In section 8 we present the method used 
to match the photometry over the whole survey in each 
band. In section 9 we show the first results on the star 
colour distributions, the galaxy colours and the galaxy 
number-counts. Einally, in section 10 we summarize the 
major steps of our photometry and we present the scien- 
tihc perspectives for the near future. 

2. OBSERVATIONS 

This survey was performed in the context of an ESO- 
key programme, thus guaranteeing the necessary observ¬ 
ing time to perform the full programme. The observations 
began in 1989 at the 3.60m telescope and were transferred 
to the NTT-3.5m in 1991. At a pace of two observing runs 
per year, the full time allocation was consumed by the fall 
1995. The observing strategy is designed to make optimal 
use of telescope time. The great advantage of the multi- 
mode instruments used (EMMI at the NTT, D’Odorico 
1990, Dekker et al. 1991, and EFOSC at the 3.6m, Buz- 
zoni et al. 1984) is the possibility to switch between spec¬ 
troscopic and photometric observations during the course 
of the same night, using the same instrument. It allows 
one to adapt to variable weather conditions : when the 
conditions are photometric with good seeing quality, pri¬ 
ority is given to imaging; if the weather degrades, one can 
switch to the spectroscopic mode by simple rotation of 
two wheels placing an aperture mask and a grism into the 
optical path. Because the number of observing nights is 
spread over many runs, the imaging data of a given run 
can be reduced in preparation for the subsequent run, thus 
yielding finding charts of the galaxies with i? < 20.5 to be 
observed spectroscopically. 

During each run, we observe our main survey region in 
first priority, and two auxiliary regions at the beginning 
and end of the night when the airmass for the main survey 
is large. Here we only report on the data for the main sur¬ 
vey region. The data in the auxiliary fields will be used as 
control samples. The position, galactic latitude and solid 
angle of the main Held are given in Table 1. In this region 
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of the sky, we acquired 25-30 frames in each of the B, V, 
R bands with the 3.6m telescope. These data cover the 
right ascension range 0^19'"00'* < R.A .2000 < 0^21’”24®, 

O O 

and correspond to an area of ~ 0.53 (i?.A.) x 0.24 {Dec). 
With the NTT telescope, we acquired 20-25 frames with 
0/121^24® < R.A 2000 < 0'‘26™00® (~ 1.00° (i?.A) x 

o o 

0.24 {Dec)). The 0.24 width in declination is covered by 
3 EFOSC Helds and 2 EMMI fields. 


Table 1. Survey characteristics 


center 


Total area 

R.A., DEC. 

(deg) 

(deg^) 

0'‘22™30b-30“06' 

-83 

0.24x1.53 


2.1. Instrumentation 

Because we performed the photometric observations over a 
period of six years, different CCD cameras were used over 
the years to follow the improvement in CCD technology, 
so that the size and the characteristics of the individual 
CCD frames vary across the mosaic (made of typically 50 
frames in each band). We summarize the principal CCD 
characteristics in Table 2. The first column gives the tele¬ 
scope and instrument used with the different channel for 
EMMI on the NTT (blue imaging channel: BIMG, and 
red imaging channel: RIMG). The useful area is the used 
part of the CCD which differs for large CCDs (LORAL, 
TEKTRONIX) from the total area because of vignetting. 

2.2. Observations 

At the 3.60m telescope with EFOSC, the photometric ob¬ 
servations were performed in a range of seeing between 
1.10 and 1.65 arcsec. When the FWHM is < 2 pixels, the 
profile of unresolved objects is poorly sampled, which de¬ 
grades the star/galaxy separation. With EFOSC, we ob¬ 
served the photometric fields with airmasses between 1.0 
and 1.4, and the exposure times for the B, V, R bands are 
respectively 30 , 25, 20 minutes. The majority of the fields 
were obtained in only 1 exposure in each band. Bright 
stars (R < 18 mag) are thus saturated. Since 1991, we 
have observed at the NTT with an average seeing of ~ 1.0 
arcsec. As the pixel size of the CCDs is smaller, we have a 
better sampling of the p.s.f. and therefore our star/galaxy 
separation is limited by seeing. The airmasses are also in 
the range 1.0 - 1.4. The exposure times for the B, V, R 
bands were respectively 25, 20, 15 minutes. To avoid satu¬ 
ration of the CCDs by bright stars, the images are summed 
from 2 or 3 exposures. The photometric strategy adopted 
to probe the different strips is to obtain a mosaic of CCD 
frames regularly offset by ^ of the CCD size providing 
many overlaps for subsequent checks of the photometry. 


In Table 2, we summarize the characteristics of the pho¬ 
tometric observations during the different runs. For each 
run, the photometric zero-points (apparent magnitude of 
a star with an absolute magnitude M = 0. and null color 
term exposed during 1 second, as defined in Equation (1)) 
are given. 

3. DATA PRE-REDUCTION 

Before the photometric analysis, the raw data are reduced 
using the MIDAS environment in the following steps: 
bias subtraction, skimming subtraction (for RCA CCD 
^8), flat-fielding and cosmic-ray removal. 

3.1. Bias subtraction 

A “master” bias frame is obtained by averaging over ~ 50 
frames from which cosmic ray events have been removed. 
As the bias level varies with time, the bias subtraction for 
each scientific frame is performed in two steps : 

(1) Subtraction of the master bias frame 

(2) Subtraction of the mean difference between the over¬ 
scan of the scientific frame and that of the master bias 
frame. 

3.2. Skimming subtraction 

With RCA CCD ^8, the raw CCD frames have a skim¬ 
ming pattern, characterized by systematic column inten¬ 
sity offsets across the entire CCD. This pattern varies 
in intensity at low illumination levels and becomes con¬ 
stant above ~ 1000 adu, well below the sky level of all 
our scientific frames but above that for the photomet¬ 
ric calibration frames. Because the skimming offsets are 
additive biases, we can calculate them by using fiat-field 
frames with different illumination levels. Scaling and sub¬ 
traction of two fiat-field frames with different exposure 
times yield a preliminary skimming frame. As the skim¬ 
ming feature is stable along each column, the signal to 
noise is improved by replacing each pixel of a column by 
the mean value along the column. The mean level of the 
entire skimming frame is then adjusted to be zero. A set 
of skimming frames derived from pairs of fiat-field frames 
with increasing illumination levels are calculated. Then 
for each science frame, the appropriate skimming frame 
can be subtracted. The skimming subtraction is impor¬ 
tant for the calibration frames where the sky level is very 
low. To preserve the quality of our photometry, the cali¬ 
bration frames for which we cannot adequately subtract 
the skimming pattern are rejected. 

3.3. Flat-fielding 

The fiat-field frame provides a map of the sensitivity vari¬ 
ations over the CCD chip. This map depends on the spec¬ 
trum of the incoming light. As the fiat-field pattern is sta¬ 
ble from night to night, we obtain for each observing run a 
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Table 2. CCD characteristics and Photometric observations 


Telescope/instrument 

GCD 

Pix size 

(7pix) 

Useful area 
(arcmin^) 

Period 

Average zero-points 

B V R 

3.60m/EFOSC 

RCA#8 

0.675 

3.7X5.7 

43-44 

23.70±0.03 

24.24±0.04 

24.21T0.02 

3.60m/EFOSC 

RCA#8 

- 

- 

45-46 

23.41±0.01 

24.06±0.03 

24.03T0.02 

NTT/EMMI (RIMG) 

THX#18 

0.44 

7.5 X 7.5 

49-50 


24.30±0.02 

24.65±0.03 

NTT/EMMI (RIMG) 

LOR#34 

0.35 

9.3X8.8 

52 

24.69±0.02^ 


25.14T0.02 

NTT/EMMI (RIMG) 

TEK#36 

0.27 

9.1x9.! 

54 


25.40T0.02 


NTT/EMMI (BIMG) 

TEK#31 

0.37 

6.2x6.2 

52 

24.69±0.02 



NTT/EMMI (BIMG) 

TEK#31 

- 

- 

54 

24.26±0.01 




“super flat-field ” in each filter band by doing a median fil¬ 
tering of all the scientific frames obtained during the run. 
The median filtering removes the objects from the images 
and yields the large scale sensitivity variations. This is the 
best flat-field frame which can be obtained, because it is 
derived from the sky on the science frames themselves. 
The number of frames must be sufficient to create a high 
signal to noise final flat-field frame. In practice, because 
our fields are sparsely populated, the super flat-held re¬ 
sults from the median filtering of 5 to 20 science frames 
(depending on the observing run). This flat-held is normal¬ 
ized and is divided into each data frame. The large-scale 
residual variations in the background of each flat-fielded 
frame are <1 %. 


3.4- Summation of multiple exposures 

After flat-fielding, we align multiple exposures of identical 
fields using several unresolved objects, sum the individual 
exposures and calculate the mean airmass for the final 
frame. 


3.5. Cosmic ray removal 

To remove cosmic rays from the summed frame, we apply 
a filtering algorithm kindly provided by P. Leisy. Each 
pixel is compared with the mean of the 5x5 neighbouring 
pixels. If the pixel value differs by more than 5a from the 
mean, it is replaced by the mean value. The value of 5cr is 
large enough to prevent the subtraction of real objects. 

4. PHOTOMETRIC ANALYSIS 

This section briefly describes the photometric software 
used in the data analysis (SExtractor : Source Extrac¬ 
tion software). A complete description is given in Bertin 
& Arnouts (1996, hereafter BA96). 

4 . 1 . Description of photometric package 
The photometric analysis is done in four steps: 


4.1.1. determination of the sky background 

The background map is derived by binning the frame into 
large meshes (32 x 32 or 64 x 64 pixels), and removing 
the possible overestimations due to bright objects using a 
K-sigma clipping algorithm. This sky background map is 
subtracted from each frame. 

4.1.2. Detection of objects 

The detection algorithm determines a group of connected 
pixels above a given threshold. A minimum number of 
connected pixels is chosen in order to avoid spurious ob¬ 
jects (bad pixels, un-removed cosmic rays,...). Convolu¬ 
tion with a gaussian with a FWHM close to the seeing 
improves the detection of objects by decreasing the back¬ 
ground noise. The typical threshold used is 1.2 ct above the 
background value ( u is measured from the unsmoothed 
background ) corresponding to a surface brightness equal 
to 27 mag.arcsec“^ in B, 26.5 mag.arcsec”^ in V and 26 
mag.arcsec"^ in R. When such a low threshold is used, 
large spurious faint objects can appear in the wings of 
objects with shallow profiles. This effect is seen around el¬ 
liptical galaxies or bright stars where the local background 
noise increases and can exceed the detection threshold. A 
cleaning procedure is therefore applied to check if the de¬ 
tected objects are real. For the faint objects in the vicinity 
of bright objects, a new estimation of the local background 
is obtained by assuming that the dominant central object 
has large gaussian wings. If the mean surface brightness of 
the faint object is lower than the calculated local thresh¬ 
old, the object is rejected. 

4.1.3. Deblending of neighbouring objects: 

Each detection of a connected set of pixels is pro¬ 
cessed through a deblending algorithm based on multi¬ 
thresholding as described in BA96. 

4.1.4. Photometry 

Three kinds of magnitudes are used: isophotal, “corrected 
isophotal” and “adaptive aperture”. 

The “corrected isophotal” retrieves the lost flux in the 
isophotal magnitude by assuming that the wings of ob- 
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jects outside the limiting isophote are nearly gaussian (see 
Maddox et al. 1990a). This estimation of the “total” mag¬ 
nitude could be improved by assuming that the profile of 
objects follows the more realistic Moffat profile. However, 
if the value of the threshold is low, the gaussian profile 
provides an adequate correction as demonstrated by the 
tests on simulated frames (see section 4.2). 

The “adaptive aperture” magnitude is the best esti¬ 
mation of the “total” magnitude. The algorithm is similar 
to the “first-moment” measure designed by Kron (1980). 
This magnitude is calculated in two steps: 

(1) The object’s light distribution above the isophotal 
threshold is used to measure an isophotal elliptical aper¬ 
ture characterized by the elongation e and position angle 
9 . 

(2) The first moment ri is calculated in an aperture 2 
X larger than the isophotal aperture in order to reach the 
light profile information below the isophotal threshold. 
The first moment is used to define the “adaptive aperture” 
of radius kri inside which we measure the “total” magni¬ 
tude. The principal axes of each object are defined by ekri 
and krije. Kron (1980) uses a circular aperture of radius 
kri with k =2 which measures 90% of the total flux from 
the objects. To converge near the “total” magnitude we 
can increase the k value, but a compromise must be found 
between the added measured flux and the increasing noise 
in larger apertures. As Metcalfe et al. (1991), we choose 
k =2.5, yielding 94% of the total flux inside the adap¬ 
tive aperture (this value was calculated using simulated 
frames with a large variety of galaxy profiles). In con¬ 
trast to the isophotal magnitude which operates at hxed 
signal-to-noise, the “adaptive aperture” magnitude may 
be determined at very low signal-to-noise. Sometimes for 
faint objects, ri may converge to erroneously small aper¬ 
tures. We therefore constrain the apertures to a minimum 
value of Rmin ~ 3.5 aab (where a at is the mean standard 
deviation of the bivariate gaussian profile defined by the 
second order moments of the object profile (BA96)). 

4-2. Simulations 

To test the quality of the photometry for objects with dif¬ 
ferent magnitudes, we use simulated frames provided by 
E. Bertin (BA96) and generated according to the follow¬ 
ing scheme. A sky background is defined by a sky surface 
brightness = 21 mag/arcsec^), and Poisson noise is 
added according to a gaussian distribution as observed 
on real CCD frames. A set of stars with a Moffat prohle 
(Moffat 1969) are generated. Galaxies are produced with 
a large variety of shapes and sizes. The pixel size and see¬ 
ing disk are adjusted to resemble as much as possible the 
real frames. The frames are defined through the R band 
corresponding to the band used here for the selection of 
the spectroscopic sample and for the star/galaxy separa¬ 
tion. The results of the tests on several simulated images 
are based on ~ 6000 galaxies and ~ 600 stars in total. 


4.2.1. Test of photometric accuracy 

In Fig. 1, we compare the mean difference between the dif¬ 
ferent measured magnitudes and the true magnitude for 
the galaxies. The mean difference is calculated in bins of 
0.5™ of the true magnitude. The error bars represent the 
r.m.s scatter around the mean. 

As expected, the isophotal magnitude which measures the 
flux inside the defined isophote looses the flux in the wings 
outside this isophote. Fig. 1 shows that an increasing frac¬ 
tion of flux is lost for fainter magnitudes. The “corrected 
isophotal” magnitude (defined in section 4.1.4) provides a 
signihcant improvement. 

Fig. 1 also confirms that the “adaptive aperture” mag¬ 
nitude (Kron magnitude) measures ~ 94% of the flux of 
objects over the entire magnitude range (the systematic 
offset of 0.™06 is indicated by the dashed line in Fig. 1). 
To reach the “total” magnitude, we then subtract a con¬ 
stant value of 0.™06 to all magnitudes. The error bars in 
Fig. 1 are ~ 0.05™ for galaxies with R < 21., and up to ~ 
0.2™ at fainter magnitudes. The interest of the “adaptive 
aperture” magnitude over the isophotal magnitude origi¬ 
nates in it insensitivity to seeing and redshift (Kron 1980). 

Because a large fraction of the detected objects in our 
survey are stars, we also compare in Fig. 2 the different 
magnitudes for the simulated stars. Even when a gaussian 
prohle is used to correct the flux lost in the wings of stel¬ 
lar objects with Moffat prohle, the “corrected isophotal” 
is very close to the total magnitude for R< 23.5. 

Because the fraction of hux measured in a hxed aperture is 
constant for stars, we check the reliability of the second or¬ 
der moments by comparing the “adaptive” magnitudes to 
the “aperture” magnitudes. The aperture radius is dehned 
by 3.5x(Ta where a a is the second order moment along the 
major axis, estimated by using unsatured stars with high 
signal-to-noise (typically 19™ < i? < 22™). The lost hux 
is stable in the full magnitude range and is close to 0.04™. 
This value is in good agreement with the “adaptive” aper¬ 
ture obtained by using a minimum radius dehned as 3.5 
aab (where aab = ^/o’a-ab is calculated for each object). 
We can then conclude that the second order moments can 
be reliably estimated down to very low signal-to-noise ra¬ 
tios and the “adaptive” magnitude is a good measure of 
the total magnitude. 

Despite the reliability of “adaptive” magnitude for 
stars, we use the “corrected isophotal” magnitude as be¬ 
ing the best magnitude to estimate the “total” magnitude 
. Because, we make the star/galaxy separation only for 
R < 22™, we use the “corrected isophotal” magnitude up 
to this cut-off for the stellar sample and the “adaptive” 
magnitude for the galaxy sample. At fainter magnitudes, 
no classihcation is done, thus we use the “adaptive” mag¬ 
nitude for all objects. 



Measured - Simulated 
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0.4 


“Kron magnitude 
0.3 p ■ Isophotal 

xlsophotal corrected 



20 21 22 23 24 

Simulated Magnitude (bin = 0.5 mag) 


Fig. 1. Mean difference in 0.5™ bins between the measured 
magnitudes and the true magnitudes for galaxies in simulated 
R frames. The symbols for the different magnitudes used are 
given inside the graph. The error bars show the r.m.s scat¬ 
ter around the mean. The dashed line represents the expected 
position of magnitude difference for the “adaptive aperture” 
magnitude corresponding to the 94 % enclosed flux. 

4.2.2. Conclusion 

The results of the tests applied to the simulated images 
show that : 

(1) For galaxies, the “adaptive aperture” magnitude 
is a robust estimate of the “total magnitude” (after cor¬ 
rection for an offset of 0.06™) down to R < 24.0. The 
uncertainty in the measure is close to 0.05™ up to R<21 
and increases to 0.2™ up to R<23.5. However, this is an 
aperture magnitude and it is sensitive to crowding by close 
neighbours. Therefore, when an object has a neighbours 
closer than 2 isophotal radii, we use by default the “cor¬ 
rected isophotal” magnitude. 

(2) For stars, the “corrected isophotal” magnitude is a 
reliable estimate of the “total magnitude” for R < 23.5. 

(3) Because we measure “total” magnitudes, the 
colours can be calculated as the difference between the 
“total” magnitudes in each passband. The colour error 
can be estimated as the quadratic sum of the error in each 
band. For objects brighter than R < 22, the colour error 
is ~ 0.07™ and at fainter magnitudes the colour error is 
~ 0.28™. 

4-3. Star-galaxy separation 

To build the spectroscopic catalogue, it is important to 
exclude the stellar objects. The star/galaxy separation 
is performed using the neural network included in the 
photometric software SExtractor. For a complete descrip¬ 
tion of the principle, the training (with using simulated 


Fig. 2. same as Fig. 1 for stars. The dashed line represents 
the flux measured inside a radius defined as 3.5X(Ta, where Oa 
is derived from stars with high signal-to-noise (as described in 
text). 


frames containing stars and galaxies), and various tests of 
this neural network, the reader can refer to BA96. Here, 
we describe the input and output parameters used for 
the separation. For each image, the classifier works in a 
ten-parameter space containing 8 isophotal areas defined 
by dividing the peak intensity in 8 levels equally spaced 
in logarithmic scale, the central peak intensity, and one 
control-parameter which defines the fuzziness of the frame 
and is chosen to be the seeing (in pixels). This set of 
input parameters is shown to provide an “optimal” de¬ 
scription of the characteristics of each image (BA96). The 
neural network provides an output parameter defined as a 
“stellarity-index”. Because the p.s.f. of a frame is the pa¬ 
rameter which determines the quality of the performance 
of the neural network, SExtractor estimates on each frame 
the FWHM of the p.s.f. using the unsatured bright stars 
and uses it for the neural network. The “stellarity-index” 
output parameter is a measure of the conhdence level in 
the classification of each image. The value varies between 1 
for stars and 0 for galaxies. In Fig. 3, we show the stellarity 
index as a function of magnitude for a R frame. For our 50 
R frames, a reliable classification can be done for R < 22 
with a success rate for the galaxies close to 95% (BA96). A 
fainter magnitude limit can be reached for several frames 
with good seeing quality (such a frame is shown in Fig. 3). 
To be sure that galaxies are not mis-classified as stars, we 
define as stars all objects with a stellarity index greater 
than 0.8. The corresponding magnitude limit for classifica¬ 
tion is 1.5 magnitude deeper than our spectroscopic limit 
of R=20.5, and thus guarantees that our spectroscopic 
sub-sample is poorly biased by mis-classified objects. 

Fig. 3 shows that at faint magnitudes R > 23, the stel- 
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larity index approaches 0.5 because no distinction can be 
made between the two classes when the profiles are dom¬ 
inated by the seeing disk. On the other hand, at bright 
magnitudes, for the largely saturated stars (R < 16) the 
stellarity index may drop below 0.8. 

When we started the spectroscopic observations, the pho¬ 
tometric software by E. Bertin was not avalaible. We 
thus used INVENTORY, the photometric software devel¬ 
oped by A. Kruszweski (West & Kruszewski (1981)) and 
avalaible in the MIDAS environment. INVENTORY pro¬ 
vides a reliable star-galaxy separation for R<21 and was 
used to generate the entry catalogue for the first half of 
the spectroscopic observations. The good success rate in 
the separation is confirmed by only 15 spectra of stars ob¬ 
served out of a total number of 521 reduced spectra (i.e. 
< 3%) (Bellanger et al. 1995a). 



14 16 IB 20 22 24 

MAGNITUDE 


Fig. 3. Star/galaxy separation for a science frame observed 
in the R band with a seeing of 0.9 arcsec. The ordinate shows 
the output parameter of neural network , the “stellarity index 
” which provides an estimate of the confidence level in the 
classification of each object as a star (CLASS = 1) or a galaxy 
(CLASS = 0). 


5. ASTRONOMICAL COORDINATES 

To obtain precise astrometry for the detected objects in 
our CCD frames, the regions of the SERC-J photographic 
plates containing our survey were digitized at the MAMA 
(“Machine Automatique a Mesurer pour I’Astronomie ” 
at Paris Observatory) which is developed and operated by 
CNRS/INSU (Institut National des Sciences de I’Univers). 

The astrometry is performed in two steps. 

(1) By alignment of the images of the digitized plates 
onto our CCD frames (using existing MIDAS commands 
applied to objects in the field), we are able to derive the 
coordinate transformation equations from CCD to photo¬ 


graphic plate. There is one set of transformation equation 
for each CCD frame. 

(2) The transformation equations from photographic 
plate to equatorial coordinates are provided by the 
MAMA facility as FORTRAN programs. The astromet¬ 
ric reduction uses the PPM catalogue (Roeser & Bastian 
1991). 

(3) Combination of the two transformations yields for 
each CCD frame the transformation equations from pixels 
to equatorial coordinates. We use a set of MIDAS proce¬ 
dures which were written with the goal to perform the 
transformation routinely for the numerous fields of the 
programme (Revenu and de Lapparent 1992). 

The overlapping regions of neighbouring CCD frames 
allow us to estimate the internal astrometric consistency. 
In Fig. 4, we plot the median position difference defined 
as ^Jcos‘^{DECo) x (RAq - RAi)^ + {DECq - DECiY 
in bin sizes of 1 mag. The NTT observations are repre¬ 
sented by dots and those from the 3m60 by crosses. The 
error bars are the 1 a dispersion measured in each magni¬ 
tude bin. For the NTT, the differences of median position 
vary from 0.1 "±0.1" for R< 22 to 0.2 "±0.3" for R<24, 
and from 0.2 " ± 0.2" for R< 20 to 0.4 " ± 0.3" for R<24 
with the 3m60. The better accuracy for the NTT fields 
is due to the better sampling (smaller pixel size) which 
improves the measurement of the central position of the 
objects: the seeing disk is spread over ~ 3 pixels for NTT 
frames, whereas it is spread over only ~ 2 pixels for 3m60 
frames. 

6. CALIBRATION 

To calibrate our CCD magnitudes into the Johnson B, V 
and Cousins R standard system (hereafter B, V, R), we 
observed the faint standard star sequence NGC 300 from 
Graham (1981) and the bright standard star sequences 
MARKA, NTPHE, PG0231, SA92 from Landolt (1992). 
Several sequences were observed each night in order to 
obtain a zero-point for each night. We use the following 
linear transformation equations to convert instrumental 
magnitudes niccd into standard magnitudes M^td '■ 

Mgtd = 'nT'ccd ~ Am X sec^ ± kM x COLOUR ± Cm 

( 1 ) 

where Mgtd define the standard magnitudes and 
niccd define the instrumental magnitudes (in adu.sec“^) 
through respectively the three filters B,V,R. The colour 
term COLOUR is the standard colour provided by Gra¬ 
ham or Landolt B—V or V—R. The instrumental CCD 
magnitude of standard stars is calculated as the “corrected 
isophotal” magnitude (see section 4.2.1), thus yielding a 
good estimation of the total magnitude. 

The extinction coefficients Am for the different filters B, 




S. Arnouts et al.: The ESO-Sculptor Faint Galaxy Redshift Survey: The Photometric Sample 



MAGNITUDE 


Fig. 4. Mean difference in arcsec in the astrometric position of 
common objects in CCD overlaps as a function of R magnitude 
in bin sizes of 1 mag. The data for the NTT and 3m60 tele¬ 
scopes are represented by dots and crosses respectively. The 
error bars show the 1 a r.m.s. dispersion. 

V, R derive from the La Silla extinction curve provided in 
the ESO manual (Schwarz & Melnick 1993). 

The zero-points Cm are specific to each combi¬ 
nation of Telescope/Instrument / CCD/Filter/atmospheric 
conditions. These are easily measured using a large 
number of standard stars observed each night (see 
Fig. 6). The colour coefficients kn allow the magni¬ 
tudes from the “observing” filter (resulting from the Tele- 
scope/Instrument/CCD/Filter combination) to be cor¬ 
rected into the B, V, R standard filters. We assume that 
they remain constant during each observing run, and use 
the calibrations of an entire run to calculate these coeffi¬ 
cients. Two color coefficients for the V band can be cal¬ 
culated depending on whether B—V or V—R colours are 
used (denoted kvs or respectively). 

In practice, we estimate the colour coefficients for each dif¬ 
ferent configuration of Telescope /Instrument /CCD /Fil¬ 
ter as shown in Fig. 5. The resulting colour coefficients are 
shown in Table 3. The quoted errors are the r.m.s. uncer¬ 
tainties in the linear fit. The measured cofour coefficients 
are then used to determine the accurate zero-point Cm 
for each night. 


7. STANDARD MAGNITUDES 

7.1. Transformation equations 

The calibration of the five different CCDs used in the 
imaging survey allows us to derive a homogeneous set of 



Fig. 5. Colour term as a function of colour for standard star 
sequences of an entire observing run (NTT/EMMI-B, CCD 
TEK#31) . The error bars for each star is the quadratic sum 
of the instrumental error resulting from several measurements 
and the intrinsic error given by the authors. The solid line 
represents the weighted least-squares regression whose slope 
provides the value of fcs. The C° are chosen for each set of 
stars of each night to minimize the scatter in the colour term 
kb. 



Fig. 6. Resulting zero-point value Cl for one observing night 
using the color coefficient ks calculated in Fig. 5. The error 
bars for stars are the same as in Fig. 5, and the zero-point is 
calculated by weighted least-squares regression by a constant 
term. 
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Table 3. Measured colour coefficients 


Inst/CCD 

ks 

kvB 


kR 

EFOSC/RCA#8 

0.16T.03 

0.04T.02 

0.10T.02 

0.00T.02 

EMMI-R/THX#18 


0.05±.01 

0.10T.02 

-O.lOT.Ol 

EMMI-R/LOR#34 




-0.03±.01 

EMMI-R/TEK#36 


0.03±.02 

0.05±.02 


EMMI-B/TEK#31 

-0.21T.02 





standard magnitudes for all detected objects. Instrumen¬ 
tal magnitudes are related to standard magnitudes by the 
transformation equations obtained from Eq.l : 

Mstd = Meed + knfiB — V) (2) 

where Meed = meed - Amscc^ + Cm 

We assume that no correction for galactic absorption 
is necessary because our fields lie at high galactic latitude 
~ —83° and the widest extent of the full survey is only 
1.53 deg. 

Because we have a priori no knowledge of the colours 
of the detected objects in the Johnson-Cousins system, we 
invert the set of equations (2) so that the colour correction 
appears in terms of observed colour. We give the equations 
for the V band where two colour coefficients are measured 


V = Veed+-, - , , .{Beed-Veed) (3) 

I - ks + kvB 

b" = Veed + , . {Veed " Rccd) (4) 

1- kvR + 

Most objects detected in one filter are identified in 
the other 2 filters (see next subsection). For the galaxies 
detected in only one filter, we derive the magnitudes in 
the other filters using the standard colours B—V =1.0 
and V—R = 0.5. These mean values are derived from the 
mean galaxy colours obtained for the survey and reported 
in section 9.2.2. To check the reliability of our colour coef¬ 
ficients listed in Table 3 and the resulting standard magni¬ 
tudes, we estimate V in two possible ways using equations 
(3) and (4). The mean differences in the 2 correspond¬ 
ing magnitudes for 7000 objects observed at the NTT and 
5000 objects at the 3m60 are < 0.005 ± 0.001 mag. This 
test demonstrates the reliability of our measured colour 
coefficients. Because the R and V frames cover almost the 
same area of the sky, the final V magnitude are calculated 
in terms of V—R colour (Equation 4). 

l.'H. Colour completeness 

Although R and V frames cover almost the same area 
of the sky, the B frames cover only about 80% of area 
observed in R and V, partly because of the smaller CCD 


size in the EMMI blue channel. Within the common area 
of the B, V and R bands, we define a colour completeness 
rate for objects selected in the B band as the fraction of B 
detected objects with associated V or R or both detections 
(in bins of 1 mag). We calculate the analogous curves for 
the V and R bands inside the same area. These various 
curves are shown in Fig. 7. Typically, the completeness 
drops below 90% 1 mag brighter than the limit of the 
catalogue in each band. At the limit of the spectroscopic 
catalogue R<20.5, the completeness in B and V is greater 
than 95%. 

8. MATCHING CCD OVERLAPS 

Because our images were taken over a 6-year period us¬ 
ing different telescopes (3.60m and NTT), instruments 
(EFOSC2 and EMMI) and CCDs ( different pixel size and 
color coefficients), it is crucial to adjust the zero-points 
of our magnitude scales over the whole survey. This is 
done by comparing the measured magnitudes of objects 
located in the overlapping edges of the CCD images (1/10 
of the CCD size, i.e. < 1 aremin). This adjustment guaran¬ 
tees an internally homogeneous photometry, and corrects 
for the systematic offsets in the photometric calibrations 
from one observing night/run to another, as well as other 
possible variations not accounted for like the systematic 
variations in the extinction curves due to the eruption of 
the Pinatubo volcano in 1991 (Burki et al. 1995). 

The input images have the standard magnitudes derived in 
section 7, and the equatorial coordinates (a, d) described 
in section 5. Because the overlaps contain few objects, we 
use all stars and galaxies with R < 22 mag. The mag¬ 
nitudes can then be directly compared in the overlaps of 
each CCD, and we show below how the mean offsets pro¬ 
vide the correction of zero-point to apply to each CCD. 

8.1. Method 

Because each CCD frame has 2 to 4 overlaps, the adjust¬ 
ment of zero-points must be done by a global least-squares 
fit. We use the method described by Maddox et al. (1990a) 
and originally proposed by Seldner et al. (1977). Seldner 
designed this method for adjustment of galaxy counts in 
the overlapping cells on neighbouring plates to correct the 
systematic variations in the magnitude limits from plate 
to plate. Instead of comparing the counts, Maddox for¬ 
malism allows to compare the differences of magnitude 
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Fig. 7. Colour completeness of B, V and R detected objects in B—V and V—R colours. The completeness is estimated in bin 
sizes of 1 mag. The coding of the curves is specified inside the graph. 


for common objects on overlapping photographic plates. 
For a complete description of the technique the reader can 
refer to Maddox et al. (1990a). Here we briefly review the 
principle of the method. 

For each galaxy k on CCD i, we define a correction for 
the zero point SCi by 

rriQ = nii + Et + SCi ( 5 ) 

where Wg is the “true magnitude” and is the mea¬ 
sured magnitude with an error On each overlapping 
CCD pair (i,j), we obtain 

rriQ = rrii + Ei + SC^ = m'^ + Ej + SCj ( 6 ) 

We define a mean offset Ty so that 

Tij =< m'; - >~ SCi - SCj ( 7 ) 

here, we assume that < Ei >=< Ej >= 0. Because 
the CCD provides a linear response in a large magnitude 
range, Ty is independent of the magnitude {rrij + mi)j2. 


We thus measure Ty by an offset in the zero*^ order of 
(mj—nii) versus (nij+mi)/2 (for the photographic plates, 
Maddox & al (1990a) define Ty by a third order adjust¬ 
ment). Fig. 8 shows an example of the differences of mag¬ 
nitudes for all objects in one overlap (i,j). We require that 
the overlaps taken into account in the equation system 
have more than 5 objects in order to guarantee a reliable 
measure of the mean offset value. 

For CCD i, we obtain several correction factors corre¬ 
sponding to each of its neighbours j which we denote as 

SCi = SCi^ = SCj + Ty j c i (8) 

where SCij is the correction factor estimated with 
neighbouring CCD j, and j G i denotes all the neighbours 
j of CCD i. 

The goal is to determine the unique value SCi for CCD 
i which globally minimizes the scatter in the SCij. These 
coefficients are calculated by minimization of the function 
F defined by: 

^ ^ E E Ti. + ^Cj - SCi)^ (9) 

i jCi 
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where Wij = l/var{Tij) are the weighting factors 
which favor the overlaps with a small variance. 

With = 0, we obtain the following set of equations : 


^Cn+l 


jdi 




( 10 ) 


which are solved by iteration. To increase the conver¬ 
gence speed and the stability of this iterative method, we 
adopt the technique from Maddox of adding the previous 
value (5(7" into equation (10): 


= 


6C^Wu + E('^^r + 

jCi 

w,, + 

jCi 


( 11 ) 


where Wu is the mean value of all Wij of each CCD i. 



CCD number 




(mj + mi)/2 


Fig. 8. Difference of magnitudes for common objects in two 
overlapping CCDs j and i. The solid line shows the linear fit for 
objects brighter than 22 mag which measures Tij = rrij — rm. 


8.2. Results 

To constrain the set of correction factors, Maddox used 
= 0. In our case, we have re-observed several fields 
across the whole survey during the last observing run: 4 
fields in B with TEK^Sl, 4 in V with TEK:;(136 and 5 
fields in R with LOR#34. The seeing conditions were good 


Fig. 9. The upper graph shows the estimated magnitude cor¬ 
rection for the R filter to be applied to each CCD zero-point. 
The first 30 CCDs are fields observed at the 3m60 and the 
remaining come from NTT observations. CCDs from 47 to 50 
are used as reference frames with no correction factors which 
provide the constraints to the other frames The lower graph 
shows the histogram of Tij before correction of the zero-point 
(empty histogram) and after correction (hashed histogram). 


{FWHM ~ 1.0", 1.0", 0.8" in B, V, R respectively) and 
the photometric calibrations were performed with great 
caution. We use these fields as fixed references for the ad¬ 
justments of the zero-points of the other frames. During 
this process, we detected systematic offsets in the initial 
V and R zero-points of the 3.6m frames derived from cal¬ 
ibration sequences (section 6) when compared with the 
NTT reference frames: SCsmeo — SCntt = -bO.lO’" in the 
R band, SC^meo — ^Cntt = —0.10’" in the V band. These 
shifts cannot be explain by a problem in the colour correc¬ 
tion because the offsets are independent of the colour of 
the objects. Because of the absence of a similar shift in the 
B band, this effect cannot be attributed to the variation 
of the extinction coefficients due to the eruption of the 
Pinatubo volcano in 1991, which occurred when the ob¬ 
servations switched from the 3m60 telescope to the NTT 
telescope. To prevent the minimization method from in¬ 
troducing a gradient in the set of correction factors around 
the boundary between NTT and 3m60 frames the equation 
systems in each band for the NTT and 3m60 are processed 
separately. 

First, we use ^ 6 Ci = 0 as a constraint for both systems. 
The solutions of the two systems provide the internal zero- 
point corrections to be applied to each CCD. Then, we use 
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the reference fields to derive the global magnitude shift of 
each ensemble of zero-points (3m60 and NTT). These var¬ 
ious steps allow us to adjust the two systems to the same 
zero-point. The resulting correction factors for the R band 
are shown in Fig. 9 a). Each point defines one CCD and 
corresponds (with some gaps) to a sequence number (see 
figure captions). The iteration of Eq. (11) is stopped when 
the difference between the previous and new estimation 
(5(7"+^ — (jC") is smaller than 0.005’” for all CCD i. The 
typical number of required iterations is 10. 

In Fig. 9 b), we plot the histogram of initial and corrected 
Tij, defined as : 

= T.,+5C^-5Ci ( 12 ) 

In these histograms, we exclude the symmetric terms 
T ■ ■ — —T 

^ ij- _ 

The r.m.s scatter in the initial Tij in B, V, R bands are 
respectively 0.11’”, 0.13’”, 0.14’”. After correction by SCi, 
we reduce the scatter to 0.04’” in all 3 bands. 

8.3. Conclusion 

These results demonstrate the efficiency of the method to 
reduce the dispersion between the CCD frames over the 
whole survey. They also allow us to estimate the dispersion 
in our magnitude system. In fact, if there was no error in 
the measures of magnitudes, Tij would be equal to zero. 
The dispersion in the after adjustments gives a 0.04’” 
error in our measured magnitudes in all 3 bands and is in 
good agreement with the 0.05’” in the adaptive aperture 
magnitude uncertainty given by the simulations (see sec¬ 
tion 4.2). To check the consistency of zero-point calibra¬ 
tions between the NTT and 3m60 fields, we have applied 
a K-S two-sample test to the B—R and B—V colour dis¬ 
tributions obtained as described in last section. The prob¬ 
ability that the 2 distributions are drawn from the same 
parent distribution is 0.45 after correction of the system¬ 
atic zero-point offsets between 3m60 and NTT frames, it 
is 0.01 when the correction is not applied. This confirms 
that correction of the offsets in zero-point for the 3m60 
in R and V bands is necessary. Finally, there could be a 
systematic shift in our overall zero-point, but it would be 
within the error bars estimated from our individual zero- 
points (i.e. < 0.05’”). 

9. FIRST RESULTS 

9.1. Star sample 

Byproducts of deep imaging surveys at different galac¬ 
tic latitudes are stellar samples in different directions of 
our galaxy at very faint magnitudes (Shanks et al. 1980, 
Kron 1980, Infante 1986, Metcalfe et al. 1991). In Fig¬ 
ures 10, II, we show for our stellar sample the different 
colour histograms B—R, B—V, for three magnitude ranges 
in V. Note that the quasars are included in this sample. 


Also, at bright magnitudes (V < 17), a large fraction of 
stars are saturated, thus raising the uncertainty in the cor¬ 
responding colours. 

At faint magnitudes {V > 18), two stellar populations 
can be seen. Brightward of U < 18, one broad blue peak 
is present near B—V~0.5-0.6. Faintward of this magni¬ 
tude a second peak appears in the red part near B—V~ 
1.4-1.5. These results agree well with the other existing 
data (see table 4 above). The red peak is interpreted as 
being nearby M stars belonging to the disk population. 
The blue sequence is interpreted as being stars belonging 
to the galactic halo (Robin et Creze 1986). 

These observations are listed in Table 4 for U <22. 
The second column gives the colours used by the au¬ 
thor, the third and fourth columns show the colours of 
the two peaks in this initial colour and the subsequent 
columns show the transformed position of the two peaks 
into the standard system Johnson-Cousins. The trans¬ 
formation equations are given in the corresponding pub¬ 
lications, and for the Kron (1980) data we added the 
transformation equations given by Majewski (1992) as 
(J—F)=0.738(B—R) -0.02. Table 4 shows that our stel¬ 
lar color distributions are in good agreement with those 
resulting from other faint CCD surveys. 

9.2. Galaxy sample 
9.2.1. Galaxy counts 
a) Counts and slopes 

We present the differential galaxy counts based on the to¬ 
tal 0.4 square degree of the survey in the R, V and B bands 
to limiting magnitudes of 23.5, 24.0 and 24.5 respectively. 
These limiting magnitudes are defined as the last magni¬ 
tude bin (Am=0.5 mag) before the turn-off in the number 
counts. To these limits the galaxy catalogs contain about 
13000, 12150 and 9500 objects in the R, V and B bands 
respectively. Note that there is no star-galaxy separation 
for R > 22 and no correction for stellar contamination is 
done because the expected number of stars at this galactic 
latitude ~ —83°)is lower than 4% (corresponding to 
an offset of ~ 0.017 of the logarithmic number counts). 

In Fig. 12, we show the differential galaxy number counts 
per square degree in 0.5 magnitude bins. The upper graphs 
show for each band the superimposed galaxy counts for all 
individual CCD fields. The dots displaced to the right of 
the columns of galaxy counts are the median counts offset 
by 0.2 mag for clarity. The error bars measure the Icr rms 
field-to-field scatter. The field-to-field scatter is ~ 20 % in 
the R band between 21 < R < 23.5, ~ 25 % in the V band 
between 21 < V < 24. and ~ 30 % in the B band between 
22 < B < 24.5. These r.m.s. dispersions significantly ex¬ 
ceed the expected Poisson variations in VN because of 
galaxy-galaxy clustering (Arnouts & de Lapparent 1996). 
In the bright part, the large values of the scatter are es¬ 
sentially due to the small number of bright galaxies per 
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B-R B-R 



B-R 


Fig. 10. B—R colour histogram for the stars selected in three different ranges of V magnitude as specified within the graphs. 




B-V B-V 



B-V 


Fig. 11. same as Fig. 10 for B—V colour. 
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Table 4. stellar colour peak for V<22 


Reference 

Observed 

colours 

blue peak red peak 

td 

1 

< 

blue peak 
B-R 

1 

> 

B-V 

red peak 
B-R 

1 

> 

Kron et al. (1980) 

J-F 

0.65 

1.75 

0.50 

0.90 

0.40 

1.50 

2.40 

0.90 

Shanks et al. (1980) 

J-R 

1.0 

2.2 

0.6 



1.6 



Infante (1986) 

J—r 

0.9 

2.3 

0.6 



1.5 



Metcalfe et al. (1991) 

{B-R)ccd 

1.0 

2.2 


1.12 



2.46 


This work 




0.55 

0.90 

0.35 

1.45 

2.40 

0.95 


CCD frame. The lower graphs give the median number 
counts. In these graphs the error bars are given as: 

^ ^ field—to—field /\/ ^field ( 1 ^) 

The differential number counts (in deg“^ 0.5 mag“^) 
are fitted by a power law in the same magnitude range as 
Metcalfe et al. (1991) for the R and B magnitudes. The 
exponent of the fit in the three bands are measured by 
a least squares fit in the logN-mag plots. These fits are 
shown by solid lines in the lower graphs of Fig. 12 and are 
parameterized as follow : 

Ngal = 10(0-460±0 02 )B ^ ^o-7.01±0.39 (^ 4 ) 

for 20.5 <B< 24.5 

Ngal = 10(0-384±0.02)y ^ io-5.03±0.47 (^ 5 ) 

for 20.0 < y < 24.0 

Ng^l = 10(O-367±O.O2)fl ^ ^o-4.42±0.47 

for 20.0 <R< 23.5 

In table 5, we summarize the results of previous CCD sur¬ 
veys on galaxy counts in the visible bands. Our slopes are 
in good agreement with the other works. 

In Fig 12, we detect two magnitude bins between 21.0 
< B < 22.0 where the density systematically decreases 
and this effect is seen in the three bands (in the inter¬ 
vals 20.5 < y < 21.5 and 20.5 < R< 21.5). First of all, 
to see if this gap was caused by inhomogeneities in the 
projected distributions, we examined the angular distri¬ 
bution (RA, Dec) in different magnitude intervals, but no 
particular feature in the clustering of the projected distri¬ 
butions was visually detected. This investigation will be 
pursued quantitatively in a forthcoming study of the two- 
point angular correlation function for these data (Arnouts 
& de Lapparent 1996). 
b) Counts comparison 

In Fig. 13, 14, 15 we compare the differential number 
counts from our data with those from the other CCD 
surveys. For the data given in other systems than the 
standard Johnson-Cousins system, we apply the differ¬ 
ent transformations provided by the authors. For the data 
from Metcalfe we apply a transformation only for the B 
band as B = Bccd + 0.15. 

For the data from Tyson, the transformation equations 


are not given. By default we use the transformations into 
the photographic system {bj.rp) given in Metcalfe et al. 
(1991) combined with the transformations from photo¬ 
graphic bands to the standard system given by Shanks 
et al. (1984). We obtain the approximate transformations 
defined as B = BTyson + 0.27 and R = Rryson — 0.07. 

For the data from Driver no transformations are done be¬ 
cause the color terms are small. Except for Metcalfe et al. 
(1991), these other works are significantly deeper than our 
data, and the number counts at very faint magnitudes (B 
> 25 and R > 23.5) are corrected for confusion. 

Figures 13, 14, 15 show that our counts in B and R are 
in good agreement with the results from Metcalfe et al. 
(1991) in both the slope and the amplitude of the loga¬ 
rithmic number counts in the red band but a small shift in 
the blue band exists as we will see below in the discussion 
of the colour distributions (section 9.2.2). The R number 
counts of our survey are significantly higher (> 3 (t) than 
those from Tyson in the common range of magnitudes. 
This difference has been interpreted by Tyson (1988) as 
being due to the a priori choice of fields devoid of bright 
galaxies, and Metcalfe et al. (1991) suggest that this differ¬ 
ence can originate from the use of isophotal magnitudes by 
Tyson (1988) in contrast to the “total” magnitudes used 
by the others authors. The data of Driver et al. (1994) also 
show a small deficit in galaxy number counts compared to 
ours at the ~ 3a level. Driver specifies that the Hitchhiker 
data suffers from a calibration uncertainty, so a small shift 
in zero-point could explain the deficit in the three visi¬ 
ble bands but this effect does not alter the slopes of the 
counts. In the data from Smail et al. (1995), the plotted 
points correspond to the average of two single fields. Their 
V counts are in very good agreement with ours but their 
R counts show a significant number excess by a factor of 
about 1.2 compared with all the others authors below R< 
24. 

Finally we compare our deep counts with the recent bright 
galaxy counts in B and R bands performed by Bertin 
& Dennefeld (1996). These counts are in good agree¬ 
ment with ours and we use them to normalize the non¬ 
evolving model kindly provided by M. Fioc. Bertin et 
Dennefeld (1996) suggest : 4)* = 0.0035hlQMpc~^mag~^ 
in the B band (estimated at B~ 19™), 4>* = 

0.0033h/^QMpc~^mag~^ in R band and we adopt an in¬ 
termediate value of 4>* = 0.0034/i5QMpc“^TOa5“^ in the 
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(a) Galaxy counts in B band 


(b) Galaxy counts in V band 


(c) Galaxy counts in R band 


Fig. 12. Galaxy number-magnitude counts by square degree per 0.5 mag interval in the B (left), V(center) and R (right) bands. 
The upper graph shows the counts for each observed CGD field and the median values with error bars given as the cr field-to-field 
fluctuation. The lower graph shows the median number count value in each bin. The error bars are given by equation (13). The 
solid line is the estimated slope by least squares fit. 


Table 5. The galaxy number count exponents for several GCD photometric surveys 


Reference 

B 

V 

R 

Tyson (1988)“ 

0.45 


0.39 

Lilly et al. (IfiSS)** 

0.38 



Metcalfe et al. (1991)“ 

0.49 


0.37 

Driver et al. (1995)“* 

0.44 

0.40 

0.37 

Small et al. (1995)“ 


0.40 

0.32 

This work® 

0.46 

0.38 

0.37 


a: in the i3j,R,I CGD photometric system 
b: in the AB photometric system 
c: in their GCD photometric system 
d: in the KPNO system 
e: in the Johnson-Cousins system 


V band {h^o is defined by Hq = 50h5okms~^Mpc~^). The 
other parameters of the luminosity function come from 
Guiderdoni and Rocca-Volmerange (1990) : a ~ —1.0 and 
M* ~ -20.6. 

9.2.2. Galaxy colours 

In this section, we present the first results of the colour 
distributions for our galaxy catalogue in the B, V and 
R bands. As shown in Fig. 7, there is a large fraction of 


objects identified in all 3 filters for B< 24 (~ 95 % of 
objects). In the following, we restrict the sample to these 
7500 common objects. 

In Fig. 16, 17, 18 we plot the mean observed colours B—R, 
B—V, V—R as a function of magnitude for the 7500-object 
sample. The solid lines draw the Icr envelope of the mea¬ 
sured colours, a being the r.m.s. dispersion of the colour 
histogram within the corresponding magnitude bin. This 
envelope is larger in B—R than in B—V and V—R because 
the expected colour in B—R varies in a large range for the 
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(a) Comparison with others deep CCD galaxy number (b) Comparison with bright photographic galaxy 
counts number counts from Bertin et Dennefeld (1996) 


Fig. 13. Comparison with others published galaxy number counts transformed into the B Johnson filter. The error bars for our 
data show the a estimated in Eq. 13. The dashed line shows the differential number counts expected for a non evolving model 
using a J?* normalized to the bright galaxy number counts from Bertin et Dennefeld (1996) as described in the text. 
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Fig. 14. Same as Fig. 13 in the V Johnson filter. 
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(a) Comparison with others deep CCD galaxy number counts 


(b) Comparison with bright photographic galaxy number 
counts from Bertin et Dennefeld (1996) 


Fig. 15. Same as Fig. 13 in the R Cousins filter 





Fig. 16. Galaxy mean colours as a function of 0.5 bin of B magnitude. The solid lines represent the Itr envelope of the measured 
colours. The upper left graph show the mean B—R from our data (filled circles) and from the data of Metcalfe et al.(1991,1995) 
(crosses). The error bars give the quadratic errors in the magnitudes obtained from the simulations. 
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Fig. 17. same as Fig. 16 for V magnitude. 



Fig. 18. same as Fig. 16 for R magnitude. The upper left graph show the mean B—R from our data (filled circles) and from 
the data of Metcalfe et al.(1991,1995) (crosses). The lower graph compare the mean V—R from our data (filled circles) with 
those from Smail et al. (1995) (crosses). 
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different galaxy types. The B—R and B—V colours clearly 
show a tendency to become bluer at fainter magnitude as 
first observed by Kron (1980) and subsequently confirmed 
by several groups (Tyson 1988, Metcalfe et al. 1991, 1995). 
At V < 21.5, the typical mean colour is B—R ~ 1.55 and 
a blueing shift to B—R ~ 1.0 is seen between 22 < V < 
24. The same tendency is visible in the B—V colour distri¬ 
bution. It decreases from B—V ~ 1.05 at V < 22 to B—V 
~ 0.6 at V = 24. Note that in Fig. 16 , at B > 24 mag, 
the completeness level drops to ~ 60% (see Fig. 7) and 
the mean colours (B—R) and (B—V) become redder due 
to the incompleteness in the R and V bands where only 
the brighter objects are identified and contribute to shift 
the colour toward redder colours. 

In contrast, the (V—R) colour in Fig. 16 - 18 shows no 
evidence of colour evolution with magnitude up to V < 
24. The same stability is obtained by Driver et al. (1994). 
In addition, we compare our observed mean colours B—R 
with those from Metcalfe et al. (1991,1995) and our mean 
V—R colours with those of Smail et al. (1995). Metcalfe et 
al’s B—R colours are systematically 0.1 mag redder than 
ours. Half of this shift is expected due to reddening in 
the fields of Metcalfe et al (Metcalfe 1995). Comparison 
with the V—R colours of Smail also show a small offset, 
but the restricted overlap in the magnitude ranges covered 
does not allow us to draw any firm conclusions about the 
agreement between the datasets. 

10. CONCLUSION 

We have obtained CCD photometry for 9500,12150,13000 
galaxies in the B, V, R (Johnson-Cousins) filters to limit¬ 
ing magnitudes of 24.5, 24.0, 23.5 respectively over an area 
of 0.4 deg^ near the southern galactic pole. Automatic pro¬ 
cedures were designed in order to reduce the data in an 
homogeneous fashion. The main steps used in obtaining 
the photometric catalogue are summarized as follows : 

i) We have applied the standard pre-reduction techniques 
of bias subtraction, fiat-fielding by median filtering and 
cosmic-ray removal. 

ii) Because different CCDs were used over the course of the 
programme, a large set of colour equations and zero-points 
were measured in order to provide reliable calibrations for 
the survey. We could then transform all instrumental mag¬ 
nitudes into the Johnson-Cousins standard system. 

iii) The survey data is a mosaic of ~ 50 CCD frames 
in each band which have ~ 1 arcmin overlaps. We could 
therefore perform an internal and global readjustment of 
the zero-point of each CCD over the whole survey. In this 
way, we reduce the zero-point field-to-field scatter from ~ 
0.11 mag to ~ 0.04 mag, in agreement with the photomet¬ 
ric uncertainty estimated using simulated images (0.05™ 
for R < 22 and 0.2™ at fainter magnitudes). We therefore 
reduce the internal dispersion in our photometry to the 
intrinsic dispersion resulting from the limitations of the 
photometric package. 


iv) Accurate astrometry is performed for each CCD frame. 
The positional accuracy is 0.1" for objects with R< 22 and 
increases to 0.3" at fainter magnitudes, as estimated from 
objects located in the overlapping edges of the CCDs. 

v) Measurement of the object colours is performed by 
matching the astronomical coordinates in each band with 
a tolerance of 1.2" in position offset. The colour com¬ 
pleteness with respect to the B band (fraction of B 
objects with R and V counterparts) is greater than 
95% up to H < 24 and drops to 60% at 24.5. The 
colour uncertainties are estimated as the correspond¬ 
ing quadratic sum of the errors in the photometric 
magnitudes {a(^B_v),(v-R),{B-R) <0.07 for R< 22 and 
^(B-v),(v-R),(B-R) <0.28 for R> 22). 

vi) The star-galaxy separation is performed using a neural 
network and is expected to be close to a 95% success rate 
for R < 22. We use the R band for the star/galaxy sepa¬ 
ration because it is our deeper band, it was observed with 
the best seeing conditions, and it is used for the spectro¬ 
scopic selection. At R> 22, no separation is done due to 
the small fraction of faint stars expected (~ 4%) at the 
high galactic latitude of the survey ~ —83°). 

The first results from the large photometric catalogue 
show that : 

i) The galaxy counts in apparent magnitude and their 
slopes in logarithmic scale are in good agreement with pre¬ 
vious CCD and photographic surveys and show an excess 
in all 3 bands with respect to the non-evolving models. 

ii) The (B—V) and (B—R) median galaxy colours show a 
blueing shift of ~0.5 mag from B ~ 22 to B ~24.5 . In 
contrast, the (V—R) median colour is nearly constant up 
to R ~ 23.5. 

iii) The galaxy counts are well fitted by a no-evolving 
model for R<20. using a 4>* derived from the new bright 
galaxy counts of Bertin & Dennefeld (1996). 

iv) Two peaks in the stellar distributions are present in 
good agreement with the Robin & Creze (1986) model of 
the galactic disk and halo. 

The next step in the study of this photometric sample 
is the analysis of the angular correlation function to {9) 
at faint magnitude (Arnouts & de Lapparent 1996). This 
work will provide information about the galaxy clustering 
up to faint magnitudes over a wide range of scales. The 
high completeness rate in colour should also allow to fur¬ 
ther characterize the properties of the faint blue galaxies. 
Together with the photometric data, the redshift survey 
of ~ 700 galaxies with R < 20.5 will provide a multi¬ 
colour optical luminosity function in the redshift range 
0.1 < 2 : < 0.5. Altogether, these data will allow us to 
address the issue of the evolution of galaxies with B < 
22. The photometric and redshift catalogues will also be 
complemented by a spectral classification of all galaxies 
with R<20.5 (Galaz & de Lapparent 1996) and will thus 
provide an unique database for studying the variations in 
the galaxy properties as a fonction of environment and 
redshift. 
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